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Abstract 

We consider a relaxation of a single mode of the quantized field in a presence of one- and two-photon 
absorption and emission processes. Exact stationary solutions of the master equation for the diagonal 
elements of the density matrix in the Fock basis are found in the case of completely saturated two-photon 
emission. If two-photon processes dominate over single-photon ones, the stationary state is a mixture of 
phase averaged even and odd coherent states. 
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In many cases, the quantum relaxation can be described in the framework of the master equation for the 
statistical operator p [1-3] (h = 1) 



dp 
dt 



H, 



IE 



2A k pA\ - A\A k p- pA{A k ^ , 



(1) 



the A k 's (k — 1, 2, • • •) being some linear operators. If the system under study is an electromagnetic field 
mode (or an equivalent harmonic oscillator), then Hamiltonian H and each operator A k can be expressed 
in terms of the annihilation and creation operators a, a 1 ^ satisfying the commutation relation [a, a^] = 1. 

There exists a specific subfamily of master equations, defined by operators A k in the form 



A k 



/i a) (n)l 1/2 a* 



1/2 



(2) 



f k a ' B \n) being arbitrary nonnegative functions of the photon number operator n = a^a. If the Hamiltonian 



is diagonal in the Fock basis, H = H(n), then Eqs 
occupation probabilities p n = (n\p\n) (n — 0, 1, . . .) 



and (0) result in a closed set of equations for the 



E [n[ +) fi a \n)p n+k 



i +) ft\n) 



n i ' '/a ■ (n-k)p n - k 



(3) 



(+) 



(n+fc)!/n!, n 



(-) 



: n\/{n— k)\. Note that Eq. (b|) does not contain off-diagonal matrix elements. 



where n 

And vice versa, the evolution of the off-diagonal elements is completely independent of the evolution of the 
diagonal ones, since the derivative d{m\p\n) /dt is expressed in terms of the elements {m+k\p\n+k) only (with 
k = 0, ±1, ±2, . . .). This means that the stationary solutions to Eqs (^)-(^) describe completely decoherent 
states, since all off-diagonal elements relaxate to zero values. 

If and flf 1 are constant positive numbers, then the terms labeled with superscripts or ' e ' describe 
the processes of fc-photon absorption or emission by some atomic reservoirs M. The choice 



fi°\n)=d k 



1 



/fc"' (n) = const 



(4) 
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corresponds to a multiphoton generalization of the Scully-Lamb |j| single-mode laser equation, the coefficient 
7fc being responsible for the saturation effect. Another equation, implying the presence of emission processes 
of all orders, described by means of powers of the shift operator up n = p n -i — p„, was proposed by Golubev 
and Sokolov ||: 

p n = r In (1 + u) p n + D[ a) [(n + l)p„+i - np n ] . (5) 

In this case, coefficients /5 e ' (n) are some rational functions determined by the Taylor expansion of function 
ln(l + w). A more general equation, with transcendental (trigonometrical) coefficients fjf (n), was obtained 

It appears that the family of known exact solutions of Eqs. (hi) and (M) is not very large. For instance, 
exact solutions for arbitrary functions /^°'(n) and fjf\ri) were found only in the stationary case with k = 1 
H . Exact time dependent solutions of Eq. Q (as well as of equations for the off-diagonal elements) in the 
case of constant coefficients and were obtained in (see also [0). Exact time evolution for the 
one-photon Scully-Lamb equation without absorption (/} = 0) was found in JxiJ - 

For fc > 2 (multiphoton processes), exact time dependent solutions of Eq. ^ with a single nonzero 
coefficient (either or fjf ) were found in Refs. [11-21]. In particular, the two-photon emission with 
= const was considered in [ p"2| , [l4[ . The case of the function (n) in the modified Scully-Lamb form 
(Q) was treated in [Jul. The two-photon absorption without emission (f^ = const) was studied in detail in 
[13—171. The case = const with an arbitrary k > 2 was investigated in ]l9| , p(i| , and the case — const 
— in | 20 [l (the evolution of the off-diagonal matrix elements in the case of two-photon absorption was studied 
in |l8f , and for fc-photon absorption — in |pp| , ^l[]). Other references can be found, e.g., in |^ , Exact 
time dependent solutions with two nonzero coefficients were obtained in [6,24-26]. In particular, the time 
dependent absorption problem with constant coefficients and was solved in Ref. [0] (a more 
detailed analysis was given recently in p?|). 

A simplified version of equation with the operator In (1 + u) replaced by the first two terms of the 
Taylor expansion, u — u 2 /2, was solved in 0, whereas the solution in a general case was given in [^5[. Other 
exact solutions with two (or more) nonzero coefficients were found in the stationary regime only. For the 
case of simultaneous fc-photon absorption and fc-photon emission (the so called systems in detailed balance) 
this was done in [B7j for the coefficients in the form and in |2i| for constant and (see also jH| 
for k — 2). A scheme of obtaining exact stationary solutions of the two-photon Scully-Lamb equation with 
single-photon losses (/{"' = const, = 0) was given in 29 1 It was generalized to an arbitrary fc > 2 in 
pb| . A stationary solution the case = An/(n + 2), — Bn(n — 1), was found in jjHJ . The case of 
three coefficients, fj^ = const, f^(n) — const or f[ e \n) = A(n + 1) _1 , was considered briefly in (3^, B3|. 
A detailed analysis of the problem with three constant coefficients, f[ a \ and (when one has one- 
and two- photon absorption, but only one-photon emission), was given recently in [B4J. 

The aim of the present article is to find a stationary exact solution to Eq. (g) in the presence of 
a two-photon emisssion. Although we did not succeed in solving the equation for a constant emission 
coefficient /j , we found that the problem can be solved in the complete saturation regime (72 ^> 1) 
of the two-photon Scully-Lamb equation (Q), when the two-photon emission is decribed by the function 
f2 B \ n ) ~ D[(n + l)(n + 2)] _1 (with the standard form — const for the absorption terms, fc = 1, 2). 
Under this restriction, there exists a 4-parameter family of equations, whose solutions can be expressed in 
terms of the confluent hypergeometric function or its special cases. 

The physical motivation for studying the new model (which is, in turn, a special case of a more general 
6-parameter family of equations admitting exact solutions) is explained by the fact that in the case of weak 
one-photon processes the stationary solutions describe an interesting class of nonclassical states, namely 
phase-averaged even and odd states (PAEOS), which are mixed analogs of the even and odd coherent (pure) 
states (EOCS) 

I \ AT l\ \ 4-1 \1 AT 2 eX P(M 2 ) M 2 ex P(H 2 ) fc\ 

\a±) = N±( a) ± - a)), N + = N_ = — (6) 

4cosh(|a|^) 4smh(|a| 2 ) 

(|a) means the Glauber coherent state 0|), introduced in |^| and studied, e.g. in [23,37-41] (for general- 
izations see, e.g. [42-46]). Since EOCS are the simplest examples of the "Schrodinger cat states" (another 
simple example is the Yurke-Stoler state [EJ \a)ys = (\ot) + i\ — a}) /s/2, the principal difference between 
EOCS and YS-states is that the EOCS have super- (even states) or sub-Poisson (odd states) photon statis- 
tics, whereas the statistics of the YS-states is exactly Poissonian), many authors considered different schemes 
of generating these states in physical processes: see, e.g. [48-52] and an extensive review [ p3| . It is known, 
in particular |5l| , that even and odd coherent states can arise due to the competition between a two-photon 
absorption and two-photon parametric processes (described by means of a nondiagonal Hamiltonian H(t)) 
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for a special initial field state. Here we show that one can obtain phase-averaged even and odd states using 
a diagonal Hamiltonian, provided that the (saturated) two-photon emission is admitted. 



2 A family of exact solutions 

A complete information about the distribution {p n } is contained in the generating function (GF) F(z,t) 
^^L Pn(t)z" . Its derivatives yield the probabilities p n and the factorial moments Af m = ^2^L m n ( n 
1) • • • (n — m + l)p„: 



Pn = -7 



1 d n F 



dz n 



Nm = 



d m F 



dz r > 



(7) 



If the products n k fl a ' (n) are polynomials of n, then one can replace the infinite system of difference 
equations (^) for the probabilities p n by a single differential equation for F(z). In the simplest cases, 
corresponding either to one-photon processes ^o[, or to a specific form of the emission operator (|^) 
fi p5|, one gets a linear differential equation of the first order. In the most of other known cases, the 
generating functions satisfy the second order differential equations of the hypergeometric type [13-18,24,26] 
(an exception is the case considered in |3l| , where a specific equation of the fourth order was solved with 
the aid of the Laplace method). One can verify that the set of stationary (p n = 0) equations (^|) results in 
the second order equation with linear coefficients {F 1 = dF/dz), 



D ( 2 a) (l + z) + [D[ a > +D 



F 



+ 



D[ a) + 2D<$ 



4 e) (l + z) + D[ e) + D$ + 3 W . 



F' 



F = 0, 



provided that functions f^'^ have the following form 



D. 



O) 



(n + l)(n + 2) : 



}[ a) (n) 



(a) 



D 



(°)» 



D[ a 2 \n + 2), 



A W (n) = £»W 



(8) 

(9) 
(10) 
(11) 



D { i j' e> and Wj° being nonnegative constant coefficients, whereas j can be any integer (excepting the value 
j — 1, which is distinguished for the sake of convenience, because it corresponds to the usual one-photon 

can be reduced to the Kummer equation |54| 

- ay = 0, (12) 



emission described by the coefficient D[ e ^). Since Eq. (|8 

xy +{c-x)y 



we have a whole family of master equations admitting exact stationary solutions in terms of the confluent 
hypergeometric function 

(a) k x k 



$(a;c;x) = ^ 



(c) k k\ ' 



(13) 



where (a)„ = a(a + 1) • • • (a + n — 1). This family is determined by 6 independent positive parameters, 
so it is larger than any one considered until now. Note, however, that we have some freedom only in the 
choice of terms responsible for the one-photon processes, while the structure of two-photon terms is fixed: 
the usual two-photon absorption and the completely saturated two-photon emission (corresponding to the 



limit 72 oo, d 2 /'y 2 -> D 2 e) in Eq. (|4j)). 

Here we confine ourselves to the special case D 



(a) 



D 



(a) 



(e) 



0. Then we have 5 independent 

parameters, an d -Dii • Normalizing all the coefficients by the two-photon absorption coefficient, D 2 a \ 

we arrive at the following set of stationary equations for the probabilities and for the generating function: 

v{(n + l)p„+i - np n - s[(n+ i)p n - np„_i] - a [p„ -p n -i]} 



-(n + l)(n + 2)p n+2 - n(n - l)p„ - r (p n -p n _ 2 ) = 0, 



(14) 
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(1 + z)F" + v{\ - sz)F' - [v(s + r 2 (l + z)] F = 0, (15) 
where we have introduced the dimensionless coefficients 

u = D\ a) /D 2 a) , s = D[ e) /D[ a \ a=D[f/D[ a \ r 2 = D 2 e) / ' D 2 a) . (16) 



A regular solution to Eq. (|15j) (without a singularity at z = — 1) satisfying the normalization condition 
F(l) = 1 can be expressed in terms of the confluent hypergeometric function 

FiA _ „m-») $("g ; ^[1 + a] ; + *]) nr . 

nj ~ *(i/ 9 ;«/[H- a ];2fl) ' (17) 

where 

R= [(^) 2 + 4r 2 ] 1/2 , ft = i(fl-i/a), ff = l[a + a + h(l + a)]. 

In particular, if s = 0, then R = 2r, h = r, and <? = jr(l + cr/r). The probabilities and factorial moments 
can be found with the aid of Eq. (Q) and the relation |54| 

— — "PI a; c; x) = — 9{a + n: c + n; x). 

dx" v 1 (c)„ v ; 

2.1 Some special cases 

If the one-photon processes dominate over two-photon ones, v —> oo (whereas r,s,a remain finite), then 
h — > 0, 7? ~ vs, ug — * 1 + a/s. Replacing (c)k by c fe in the Kummer series ( |T3| ) for c > 1, we obtain the 
following limit of the generating function (y_7|) as i> — > oo: 

/ 1 — s \ 1+<T / S 

= \T TZl ■ ( 18 ) 



This is the GF of the negative binomial distribution, which was considered in connection with the problems 
of quantum optics, e.g. in Refs. |5^, |H^]. For a — (|l8|) becomes the GF of the thermal (Planck's) 
distribution, whereas for s — > it goes to the GF of the Poisson distribution. Evidently, Eq. (hs|) is valid 
only for s < 1, whereas the general formula (jnj) holds for all nonnegative parameters s,a,u,r. If s > 1, the 
asymptotic behaviour of GF at v 1 is more complicated. For instance, in the particular case r = a — 0, 
the distribution {p n } becomes Gaussian when v 3> 1 and s > 1 J34| . 

The Poisson distribution arises also in the limit s — + oo. Then the GF tends to exp(,z — 1), i.e. the limit 
distribution has the mean photon number n = 1, independently on the values of other (finite) parameters, 
v, a,r. 

Another simple expression for the GF is obtained in absence of the two-photon absorption, D% = 0. 
Then instead of Eq. (|l^) we get the first order equation 

u(l - sz)F' - [u(s + a)+ p(l + z)]F — (19) 

(p = D 1 ^ / D^), whose normalized solution reads 

F(z)= (I^iL) 1+7 exp^(l- 2 )l , 7 = -(* + />+-)• (20) 
VI — szj is ] s \ s J 

3 Phase-averaged even and odd states 

Now let us consider the situation, when the two-photon processes dominate over the one-photon counterparts. 
Suppose first that the one-photon processes are completely absent, i.e. v — vs = va — 0. Then the solution 
of Eq. ( p"5| ) satisfying the condition -F(l) = 1 reads 

F(,) = (l-/3)^M+/3^, (21) 
cosh(r) sinh(r) 

so the occupation probabilities are given by 

_ (l-/3)r 2fc _ (3r 2k +^ 

P2k ~ (2fc)l cosh(r) ' P2k+1 ~ (2k + 1)! sinh(r) ' K ' 

The distribution ( p^ ) is nothing but a combination of the photon distribution functions of the even and 
odd coherent states \a+) and \a~) (|(J) with relative weights 1 — f3 and /3, respectively, provided that \a\ 2 



4 



is identified with the ratio of the two-photon emission and absorption coefficients r. The relative weight of 
the odd states /3 is determined by the initial conditions, since there is no correlation between even and odd 
states: /? = Sj*L P2fc+i(0). 

Using the known Wigner function of the Fock state |n)(n| [ p7| , p8| 

W n (p, q) = 2(-l) n exp (V - q 2 ) L n (2p 2 + 2q 2 ) 

(L n (x) being the Laguerre polynomial) and the generating function of the Laguerre polynomials |^|] 

^ — L n (x) = e z Jo (2^/xz) , 

it is not difficult to write an explicit expression for the Wigner function of the mixed state p = y"]p n |n)(n| 
with the coefficients given by Eq. (|2^): 

W(p,q;(3,r) = ^j'^r) " { t 1 ~ (1 ~ ^"'l J ° (^ 8r (P ' + g2) ) 

+ [(1 - 2/3)e 2 '- - f] Jo ( VM^+g^) } ■ (23) 

Here Jo (z) is the Bessel function and Io(z) is the modified Bessel function. The Wigner function ( p^ ) has 
zero mean values of the quadratures q and p, and it is very different from the Wigner functions of the pure 
even/odd coherent states 

W±(p,q;p,q) = 2iV± j exp [- (q - q) 2 - (p - p) 2 ] + exp [- (q + q) 2 - (p + pf] 

± 2exp(-q 2 -p 2 )cos[2(gp-pg)]|, (24) 

where p, q are the mean values of the quadratures in the initial coherent state \a) with a = (q + ip) /y/2. 
However, assuming q = v^cos^, p — \/2rsini/j and averaging W±(p, q;p, q) over the angle ip according to 
the formula 

, r2n , 

W±(p, q;r)= I — W±(p, q; V2rcos<p, V2r sin ip) 
Jo 2n 

we arrive exactly at Eq. (^) with /3 = for the even states and /3 = 1 for the odd states. For this reason 
we call the state described by the Wigner function (|23|) as the phase-averaged even/odd state (PAEOS). The 
phase- averaged coherent states were considered in [ |59| in connection with the problem of a classical limit for 
the quantum oscillator. The PAEOS are quantum mixtures, since the purity coefficient 



/' 



= T ^ - \ { (c^) 2 ^ + J ^ + GmL) 2 ^ ~ J ^ 



is less than 1 for r > 0. It is a monotonous function of r, whose asymptotics are 

M ^ (I-/?) 2 (I-r 2 ) + /3 2 (I-r 2 /3), r « I, (i « [(I - /3) 2 + {3 2 } r » 1. 

Nonetheless, PAEOS are nonclassical states, since the Wigner function W(g,p) (p^) assumes negative values, 
as shown in Figs 1 and 2, where we plot W(q,p) as a function of x = y/q 2 + p 2 . If r > 1, then the 
plots c orres ponding to parameters 1 — /? and /3 have a mirror symmetry with respect to the x-axis for 
x < yJr/2, since in this region the contribution of the oscillating function e r Jo (xy/8r) is dominating 
(note that W(0,0;r, (3) — 2(1 — 2(3) does not depend on r). However, the dependence on /3 disappears 
for x > y^r/2 > I, where the Wigner functions are close to zero in a wide interval (which increases with 
an increase of r), then exhibit wide and not very high maxima (at x w \/2r), and finally tend to zero 
exponentially for x 3> \/2r. In the special case /3 = 1/2 the Wigner function (G3J) is positive and does not 
oscillate even for large values of the parameter r, as shown in Fig. 3. Note that the purity coefficient fi also 
attains its minimum (for a fixed value of r) when /3 is close to 1/2. 

The type of the photon statistics (sub- or super-Poissonian) is determined by Mandel's parameter Q = 
A/2/A/1 — A/i. In the case of PAEOS this parameter equals 

Q = ^ (l - B 2 ) , B = (1 - P) tanh(r) +/3coth(r), 
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so the photon statistics becomes sub-Poissonian for /3 > i (1 — e~ 2r ), i.e. 1 — 2/3 < e~ 2r . In particular, in 
the case /3 — 1/2 (Fig. 3) we still have the sub-Poisson statistics. 

Till now we assumed that we had no one-photon processes at all. Now let us allow a small (but nonzero) 
coefficient v (weak one-photon processes). Then we can simplify Eq. (1171) with the aid of the relation 
lim„_,o <$>(av\ cv\ x) = 1 + (a/c) (e x — 1). In this limit, R — »• 2r, h — »■ r, g — > 1^(1 + s) + (s + a)/(2r), we arrive 
again at Eq. (Ell). The essential difference is that now the coefficient /3 is determined not by the initial 
conditions, but by the relative strengths of the emission and absorption processes: 

_ sinh(r) sinh(r) + (S/r) cosh(r) sj-a 
- suih (r) cosh(2r) + (S/r) sinh(2r) . * - s + 1 . I 25 ) 

and it is always less than 1/2, since 

1 - 2/3 = [cosh(2r) + (S/r) sinh(2r)] _1 > 0. 

It is remarkable that parameter v does not enter the formulas describing the stationary distribution in the 
limit v <C 1. It influences only the relaxation time t re i ~ v 1 , but not the form of the stationary state. 
If r (no two-photon emission, two-photon absorption only), then f3 — > 5/(1 + 25). If we have no 
one-photon emission (5 — * 0), then (3 — > tanh 2 (r)/ (l + tanh 2 (r)). The maximal value /3 = | is achieved 
when 5 — > oo or r — > oo. 

Mandel's parameter reads now 

= r[l-(5/r) 2 ] [l-tanh 2 (2r)] 

^ [1 + (S/r) tanh(2r)] [(S/r) + tanh(2r)] ' { ' 

The photon statistics is sub-Poissonian if r < 5, and super-Poissonian if r > 5. For a fixed r, function Q(5) 
monotonously decreases from Q(0) = rcoth(2r) [l — tanh 2 (2r)] to Q(oo) = — rcoth(2r) [l — tanh 2 (2r)] . 
Consequently, — | < Q(r,S) < |. 



4 Conclusion 

Let us formulate the main results of the paper. We have found an exactly solvable 6-parameter family of 
stationary master equations for the diagonal elements of the 1-mode field in a cavity in the presence of 
competing one- and two-photon emission and absorption processes, and we gave explicit solutions for its 
4-parameter subfamily describing the completely saturated two-photon emission regime. We have shown 
that in the limit case of weak one-photon processes, the field mode goes to the nonclassical stationary state 
which can be considered as a mixed analog of even and odd pure coherent states. Although we considered an 
idealized case of a completely saturated two-photon emission, the results obtained could help to understand 
the qualitative features of real (partially saturated) processes. 
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Figure 1. Wigner function W(x), x = \J q 2 + p 2 , of the phase-averaged even state {{3 = 0) for r = 10. 



Figure 2. Wigner function W(x), x = y/ 'q 2 + p 2 , of the phase-averaged odd state (j3 = 1) for r = 10. 



Figure 3. Wigner function W(x), x = yj q 2 +p 2 , of the phase-averaged "maximally mixed state" 
(/3 = 0.5) for r = 10. 
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